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I. INTRODUCTION 



Supersymmetry (SUSY) provides an interesting framework for stabilizing the electroweak 
scale against radiative corrections. A realistic realization of this idea is the minimal super- 
symmetric standard model (MSSM), the model with minimal field content in which all the 
standard model particles have their superpartners. In the MSSM, however, the mecha- 
nism for SUSY breaking is not specified, but is simply parametrized by introducing a large 
number of soft SUSY breaking parameters by hand. The MSSM is therefore best viewed 
as the low energy effective theory of some more fundamental theory, which must give a 
phenomenologically viable explanation of the origin of SUSY breaking. 

Models with gauge-mediated SUSY breaking |l|J^ are candidates for such a theory, and 
provide an appealing way to generate SUSY breaking soft terms in the MSSM Lagrangian. 
In these models, SUSY breaking originates in a dynamical SUSY breaking sector and is 
then mediated via gauge interactions to the MSSM sector. The SUSY breaking mechanism 
does not distinguish between flavors, and so, in the MSSM sector, all scalars with identical 
gauge quantum numbers have the same SUSY breaking mass (at the scale where the SUSY 
breaking fields are integrated out). Thus, the universality of squark and slepton masses 
at some scale is guaranteed, and serious SUSY flavor-changing neutral current (FCNC) 
problems can be evaded. 

To preserve the natural suppression of FCNCs in gauge-mediated SUSY breaking models, 
it is essential that no large intergenerational scalar mixing terms be introduced by super- 
gravity contributions. In the absence of assumptions concerning the flavor structure of such 
supergravity contributions, these mixing terms are expected to be of order the gravitino 

1 /2 

mass m3/2 ~ Fdsb/Mp, where -Fdsb is the scale of dynamical supersymmetry breaking, and 
Mp is the Planck mass. FCNC constraints then require ^3/2 <^ my/. The gravitino is 
therefore the lightest supersymmetric particle,^ which has a number of important implica- 
tions for collider phenomenology. First, the lightest standard model superpartner is now the 
next lightest supersymmetric particle (NLSP) and eventually decays into its standard model 
partner and the gravitino. The NLSP may then be charged or colored, as astrophysical and 
other constraints on such particles are thereby evaded. Second, the NLSP may be stable or 
unstable in collider experiments, depending on Fdsb- If -^dsb is high enough, decays to the 
gravitino is suppressed, and the NLSP decays outside the detector. 

The possible collider phenomena in gauge-mediated models are largely determined by the 
character of the NLSP. Typically, in these models the NLSP is either the lightest neutralino 
or the stau mass eigenstate As this particle may be either long- or short-lived there are 
four broad possibilities. If the NLSP is a long-lived neutralino, the SUSY signatures are the 



We assume that there are no other exotic light particles. 

^In the latter case, the role of the NLSP may be shared by the sleptons e-R and fiR if they are 
sufficiently degenerate with fi. 

^Throughout this paper, "long-lived" and "short-lived" refer to particles that typically decay 
outside and inside collider experiments, respectively. 
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conventional I^t signals and are identical to those present in gravity-mediated models with 
/2-parity conservation. If the NLSP is a short-lived neutralino, superparticle production 
typically results in hard photons, a signal that has also been studied in great detail 
On the other hand, the NLSP could be the stau. In the case of a short-lived stau NLSP, 
there are a number of interesting new signatures, which have been studied in Refs. 

Here, we explore the remaining possible scenario, in which the NLSP is a long-lived stau. 
This scenario is in fact realized in large portions of the parameter space: as we will discuss 
in Sec. ||, the stau is often the NLSP in models with non-minimal messenger sectors, and 
it is stable for -F^gg > 10^ GeV, as is also the case in many models, including, for example, 
the gauge mediation models of Ref. [Q. In this scenario, the familiar and hard photon 
signals are absent. However, SUSY may manifest itself at colliders through a number of 
spectacular new phenomena. For example, staus, as massive charged particles, may be slow 
and appear as highly ionizing tracks 0,§J^J^. If, on the other hand, the staus are quite 
relativistic, they cannot be distinguished from muons, and we may have excesses of dimuon 
or multi-lepton events as a result of superparticle production. 

In this paper, we present a detailed discussion of the prospects for SUSY discovery in 
this scenario, and estimate the discovery reach for a variety of signatures. At present, the 
most stringent bound on these scenarios comes from searches for stable charged particles at 
LEP [^]. Such searches, using combined data from runs with center of mass energies up to 
^/s = 172 GeV, yield the limit rrif^ > 75 GeV |IlO|, where no degeneracy between generations 
is assumed. In this study, we begin by considering possible probes from Tevatron Run I data 
with ^/s = 1.8 TeV and an integrated luminosity of L = 110 pb~^. Working within a specific, 
well-motivated model, we find that this currently available data probes regions of parameter 
space beyond the current LEP bound and even beyond the final LEP II kinematic limit. We 
then consider possible improvements at Tevatron Run II and the possible TeV33 upgrade. 
As a number of signals are essentially background-free, the prospects for improvement are 
especially bright. For ^/s = 2 TeV, and integrated luminosities of L = 2, 10, and 30 fb~^, 
we find that SUSY may be discovered for right-handed slepton masses of 110, 180, and 
230 GeV and Wino masses of 310, 370, and 420 GeV, respectively. While exact discovery 
reaches must await detailed experimental analyses beyond the scope of this study, these 
results indicate that the prospects for such searches are promising, as future Tevatron runs 
will extend current search boundaries far into ranges typically expected for superpartner 
masses. 

The organization of this paper is as follows. In Sec. 0, we briefly review the framework 
of gauge-mediated supersymmetric models. We describe a particular, well-motivated model 
that displays many of the generic features of gauge mediation models, and which we will 
use throughout the rest of the paper. In the subsequent sections, we explore three new 
signatures for long-lived staus. For each, we first estimate the discovery reach in terms of 
physical slepton and gaugino masses. In Sec. |IT1|, we discuss the possibility of detecting the 
highly ionizing tracks of heavy charged NLSPs. In Sec. |I^, we consider a possible excess of 
dimuon-like events due to the pair production of fi . The dimuon excess is a relatively weak 
probe at the Tevatron; we discuss also its prospects at the LHC. We then discuss the multi- 
lepton signal in Sec. |V|. Finally, in Sec. [V^ we summarize all of these results by studying 
and comparing the discovery reaches of the various signals in the context of the model we 
consider and its fundamental parameter space. We present our conclusions in Sec. |VI1| . 
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II. GAUGE-MEDIATED MODELS AND THE LONG-LIVED STAU NLSP 



The probes discussed in the following sections have relevance for many models, both 
supersymmetric and non-supersymmetric, with charged particles that are stable within the 
detector, and thus they have wide applicability. Even within the framework of SUSY, such 
particles are possible in a variety of settings. However, among SUSY scenarios with long-lived 
charged particles, by far the most natural are models with gauge-mediated SUSY breaking, 
where the NLSP may be charged and decays to the gravitino with macroscopic decay lengths 
for a large range of parameters.^ In this section, we will discuss these gauge-mediated SUSY 
breaking scenarios, and determine what parameter ranges lead to the long-lived charged 
particle scenario. 

In the general framework of gauge mediation, any of the MSSM superpartners could 
in principle be the NLSP and stable in collider detectors. For concreteness, and to obtain 
quantitative results, we will specify a particular, well-motivated model. In this section, we 
briefly review this simple and calculable model, which we will use to determine the promise of 
various charged NLSP signatures. In doing so, we define our conventions and notation, and 
introduce the parameters that will play a central role in the following discussions. Further 
details of this model may be found in Refs. PJTT . 



In the class of gauge mediation models we consider, there are three sectors: the dynamical 
SUSY breaking (DSB) sector, the messenger sector, and the MSSM sector. In the first, the 
F component of some chiral superfield condenses, generating a vacuum expectation value 
-^DSB 7^ 0. Fdsb determines both the mass of the gravitino 

-^DSB 

where = Mp/^/Sir ~ 2.4 x 10^^ GeV is the reduced Planck mass, and the strength of its 
interactions through, for example, 

2 

T7l~ — 

Ant = T7^^W;i + h.C. , (2) 
-TDSB 

where is the longitudinal component of the gravitino, the Goldstino. 

This dynamical SUSY breaking is mediated by two-loop effects to the messenger sector, 
which contains a singlet field S, and, we will assume, pairs of vector-like 5-1-5 rep- 
resentations of SU(5), that is, vector-like representations of SU(2)i doublets, and 
vector-like representations of SU(3)c triplets, where SU(2)i and SU(3)(7 are the standard 
model gauge groups. (Such messenger field content preserves gauge coupling unification. 
Vector-like representations 10 -|- 10 of SU(5) may also be present, and each is equivalent to 
3 vector-like 5-1-5 representations for the following discussion.) The scalar and auxiliary 
components of the singlet field then get vacuum expectation values (S) and Fs, respectively, 



''Long-lived charged particles are also possible in gravity-mediated models if the decays of a 
charged superparticle, such as a chargino or charged slepton, are, for example, highly phase space 
suppressed, or are possible only through small i?-parity violating couplings. 
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which generate SUSY breaking masses for the vector-hke fields. We denote the scale of 
these masses, the messenger scale, by M. 

Finally, once the messenger fields are integrated out at the scale M, soft SUSY breaking 
masses for the superparticles in the MSSM sector are generated. Letting 

A-^, (3) 

and assuming A <C M, we find that one-loop diagrams induce gaugino masses 

M,{M)=N,Ac,^^, (4) 

where i = 1,2,3 for the U(l)y, SU(2)l, and SU(3)c groups, ci = |, and C2 = C3 = 1. 
Two-loop diagrams induce the soft SUSY breaking sfermion squared masses 



3 r 9/,^Nn2 



i=l 



9KM) 



(5) 



Here, C( = with Y = Q — being the usual hypercharge, and C- = for gauge 
singlets, | for SU(2)i doublets, and | for SU(3)c triplets. Finally, trilinear scalar couplings 
{A terms) are also induced, but are highly suppressed, as they are generated at two-loops, 
but have dimensions of mass, not mass squared. They may be taken to vanish at the scale 
M. 

Once the boundary conditions are given at the messenger scale M, SUSY breaking pa- 
rameters at the electroweak scale can be obtained by renormalization group (RG) evolution.^ 
In particular, the Higgs mass squared can be driven negative by the large top Yukawa cou- 
pling constant, and hence break electroweak symmetry. In our analysis, we do not specify 
the mechanism by which the supersymmetric Higgs mass /z and SUSY breaking Higgs mix- 
ing mass m| are generated. Several attempts to explain the origin of these parameters may 
be found in the literature Here we regard them as free parameters. We then fix 

one combination of them by demanding the correct value of the Fermi constant using the 
tree-level Higgs potential. The remaining freedom may be specified by choosing tan j3, the 
ratio of the vacuum expectation values of the two neutral Higgs scalars, and sign(/i), where 
we follow the conventions of Ref. [0 in defining the sign of /i. Thus, the simple model we 



have defined has 4-1-1 free parameters:^ 

A^5 , M , A , tan /? , sign(/i) . (6) 



^Throughout this study, we use one-loop RG equations unless otherwise noted. 

^If all the couplings in the model are 0(1), (S) and F^^'^ (and hence M and A) are of the same 
order. However, in general, there may be a hierarchy between them. For example, if the coupling 
constant for the term in the superpotential is small, and all other couplings in the superpotential 
are 0(1), {S) is enhanced relative to Fg, resulting in A ^ M [^. Thus, in our analysis, we treat 
M and A as independent parameters. 
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In this model with the radiative electroweak symmetry breaking condition, the mass eigen- 
values and mixing parameters for all superparticles are determined once we fix these param- 
eters. 

While the determination of SUSY parameters at the electroweak scale requires detailed 
calculation, it is useful to discuss the qualitative features of the superparticle spectrum. In 
the model described above, usually /i > Mi,M2. Furthermore, as can be seen in Eq. (|[), 
gaugino masses are proportional to the corresponding gauge coupling constants, and hence 
Ml < M2. As a result, Xi ~ -B with mass m^o ^ Mi, and X2 ~ Xi ~ with 

masses m^o,m + ~ M2. In the scalar sector, sfermion masses are primarily proportional to 
the relevant gauge coupling constants. Therefore, squarks are typically much heavier than 
sleptons, and right-handed sleptons, which have only U(l)y quantum numbers, are the light- 
est sfermions. The sleptons cr and jln are almost degenerate, as Yukawa coupling effects 
are negligibly small. On the contrary, the tau Yukawa coupling may significantly lower the 
lighter stau mass through both RG evolution and left-right mixing. These effects are en- 
hanced for large tan /3. The lightest scalar is thus the lighter stau fi, which is predominantly 

We are interested in the case where < m^o and fi is the NLSP. This relation is 
most strongly and obviously dependent on N^. As one can see in Eqs. (^) and (^, at the 
messenger scale, Mj cx: N^, while Mj oc -y/iVs; thus, large reduces mf^/m^o. There is 
also an important dependence on M, as RG evolution increases right-handed slepton masses 
relative to gaugino masses, and so large values of M increase m^j/m^o. As noted above, 
large values of tan/3 decrease m^^ and so reduce mfj/m^o. Finally, the ratio mfj/m^o is 
independent of A and sign(/i) to the extent that the relation mf^/m^o ^ Mf^/Mi is valid, 
where Mf^ and Mi are the soft SUSY breaking parameters. 

In Fig. [I|, we show contours of Mf^(MMssM) = ^i(^mssm) in the (M, A^s) plane for 
various tan/3 and Mmssm = 1 TeV. (See also Ref. fl^l-) The shaded region is excluded by 
the requirement that the gauge coupling constants remain perturbative up to the GUT scale 
Mqut — 2.0 X 10^^ GeV under two-loop RG evolution. It is important to note that M; 



and Ml are only approximately the physical fi and Xi masses because of mixing effects and 
D-term contributions, and, in particular, the difference between Mf^ and m^^ may be large 
for high tan (3. However, this figure gives a qualitative picture of the parameter region where 
fi is the NLSP. For a minimal messenger sector with A^s = 1, fi is not the NLSP (unless 
tan/3 is large PJTT[]). However, for > 2, fi is the NLSP for a wide range of parameter 
space. 

In Figs. and |], we give contour plots of the soft SUSY breaking parameter ratios 
Mf^/Mi and Mj^/M2, respectively, where these ratios are at the scale Mmssm = 1 TeV. 
We fix tan/3 = 3 for these plots. In the stau NLSP region allowed by perturbativity, 
Mf^/Mi > 0.6, and for a significant portion of this region, the NLSP and Xi ^^^e fairly 
degenerate. We see also that, in this region, it is almost always true that Mj^ < M2. In 
addition, M|-^/M2 > 0.7, and for a significant portion of this region, the II and Wino states 
are also fairly degenerate. These observations will be of use in the following sections. 

Given the interaction of Eq. (El), the fi decay width is 
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and its decay length is 



L ~ 10 km X (/37) 



DSB 



10^ GeV 



100 GeV 



(8) 



where p is the n velocity and 7^(1- (3"^)-^/^. For F^^^b ~ 10^ GeV, the NLSP may be 
considered stable in detectors, as the likelihood of an NLSP decaying within the detector is 
negligible for Tevatron luminosities and the typical allowed ranges of supersymmetric cross 
sections. Such values of -F^gB are expected in many proposed models. For example, in 
typical messenger sector models, -Fdsb is likely to be larger than ~ 10'' GeV if there is no 
fine-tuning; otherwise, the SUSY breaking scale in the MSSM sector is so low that slepton 
masses become lighter than the experimental limit [l^ . There are also many other recently 
proposed models, including those based on direct gauge mediation, that give high SUSY 
breaking scales [p!4| , |T6| .Pl In a wide variety of models, therefore, the stau decays outside the 
detector, and we will concentrate on this possibility for the rest of this study. 



III. HIGHLY IONIZING TRACKS 

Long-lived staus are charged, weakly interacting, penetrating particles, and so appear in 
the tracking and muon chambers of collider detectors with little associated energy deposition 
in calorimeters. They therefore look like muons. However, since long-lived staus must have 
mass above 75 GeV |^, their velocity is typically not large if produced at the Tevatron. 
Moderately relativistic charged particles lose energy in matter primarily through ionization. 
The rate of energy loss —dE/dx is determined largely by the particle's velocity /3, grows 
rapidly with decreasing (3 for low (3 |18], and may be measured in tracking chambers. As a 
result, the cross section of slow staus may be significant, and highly ionizing tracks may be 
a spectacular signal of such particles. 

If the stau is the NLSP, all superparticles eventually decay to staus, and so all su- 
perpartner production processes may result in highly ionizing tracks. However, in typical 
gauge-mediated scenarios, superparticle production at the Tevatron is dominated by only 
a few processes. First, squarks, Higgsinos, and heavy Higgs bosons (which may decay to 
superparticles) are too heavy to be produced in significant numbers. Second, left-handed 
sleptons, though typically much lighter than these particles, are still significantly heavier 
than right-handed sleptons, and so, in terms of overall rate, the Drell-Yan production of 
left-handed sleptons is negligible relative to that for right-handed sleptons. The dominant 
processes are therefore right-handed slepton production and gaugino production. In the 
following subsections, we discuss these two processes and their potential for SUSY discovery 
in turn. 



'''Note, however, that in a class of models without messenger sectors [17], -F^gg can be as small as 
10^ to 10^ GeV, and the NLSP may decay inside the detector. 
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A. Slepton Production 



At the Tevatron, direct right-handed slepton production takes place through the Drell- 
Yan processes 

pp->l\Z* ^ fif*, crc*^, jxRix*^ . (9) 

Cross sections for a single generation of right-handed slepton production as a function of 
slepton mass are given in Fig. ^ for ^/s = 1.8 TeV and = 2 TeV.0 Here and in the 



following, we use MRS Rl parton distribution functions pO[ and the BASES phase space 



integration package [^. The cross sections for the different generations are nearly identical 
for small tan /3, since then fi ~ tr and the three sleptons are highly degenerate. In the cases 
of cr and fiR pair production, the directly produced slepton may either be stable within the 
detector or may decay through virtual Binos via Ir Irtti if kinematically allowed |^.| 
However, in all cases, the signal is a highly ionizing track from a long-lived charged slepton 
with no hadronic activity within some isolation cone.|3 

Highly ionizing tracks may be faked by a real muon that overlaps with another particle 
track, which enhances the measured —dE/dx. The role of the muon may also be played 
by a hadron that punches through to the muon chambers. However, such backgrounds 
typically arise when the muon or hadron is part of a jet, and so are drastically reduced 
by requiring the highly ionizing particle to be isolated from significant hadronic activity. 



CDF has reported such a search for highly ionizing tracks. They require the slow 



particle track to have \ri\ < 0.6 so that it is detected in the highly shielded central muon 
upgrade detector. It is also required to have —dE/dx corresponding to 0.4 < /37 < 0.85. 
The lower limit ensures that the stau reaches the central muon detector. The upper limit 
corresponds roughly to the requirement that the track be at least twice minimally ionizing. 
Such cuts are similar to those imposed in LEP analyses to eliminate background 0. They 



are also expected to reduce the background at the Tevatron to neghgible levels , and are 
appropriate for future Tevatron upgrades [^. We will assume that after all of these cuts, 
the signal is essentially background-free.|^ 



^Next-to-leading order effects fiave been shown recently to enhance the tree-level cross sections 



by 35% to 40% [19|, and so improve our discovery reaches slightly. 



^For a rather narrow range of parameters with nif-^ < m^o < ""T-f^, decays Ir IrXi ^'Iso 
possible [^], but we will not consider this scenario here. 

^°If the Ir decays in the detector with a macroscopic decay length, one might hope to find evidence 
for Ir production by searching for leptons with displaced vertices or for kinks in tracks formed by 
a Ir and its fi decay product. Unfortunately, the requirement that the decay be sufficiently far 
from the interaction point typically requires Ir-ti mass splittings so small that the leptonic decay 
products are too soft and the kinks are too small to be detected. See, for example, Ref. M. 
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Note that measurements of = \p\/m from —dE/dx and pT from the tracking chambers 
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We may now estimate the number of slepton pair events resulting in highly ionizing 
tracks. In Fig. ^ we plot the cross section for fi pair production for ^/s = 2 TeV, where 
we have required at least one fi with \rj\ < 0.6 and /37 < 0.4, 0.7, 0.85, and oo. Results 
for a/s = 1.8 TeV are similar, and we omit the corresponding figure. We find that the 
requirement (3^ > 0.4 has little effect on the signal. The efficiencies for satisfying the upper 
bound criterion < 0.85 are 25%, 44%, and 65% for m^, = 100, 200, and 300 GeV, 
respectively. While this degrades the signal substantially, slow stau cross sections on the 
order of 0.1 fb are still possible at slepton masses of 200 GeV. 

Finally, we determine the discovery reach for long-lived charged sleptons. We sum over 
three generations.^ To crudely simulate the effects of geometric triggering efficiencies and 
track quality requirements, we have included an assumed experimental detection efficiency of 
75%. With the assumption that the cuts described above reduce the background to negligible 
levels, we require, for a given luminosity, a cross section after including cuts and efficiencies 
corresponding to 5 or more events for discovery. The resulting reaches for representative 
Tevatron data samples are 

nii^ = 50 GeV (v^ = 1.8 TeV, L = 110 pb^^) 
mj = 110 GeV (v^ = 2 TeV, L = 2 fb"^) 



m 



Ir 



180 GeV (v^ = 2 TeV, L = 10 fb"^) 



m 



= 230 GeV (v^ = 2 TeV, L = 30 fb"^) . (10) 



We see that the discovery reach of 50 GeV for the currently available Run I data is not 
competitive with current LEP analyses. However, for future Tevatron runs, the discovery 
reach will be competitive with and eventually surpass the kinematic limit of LEP II. 



B. Gaugino Production 

Aside from right-handed slepton production, the other major SUSY production mecha- 
nisms at the Tevatron areQ 

pp-^W* xfxl and pp ^ 7*, Z* xfxi , (H) 

where the neutralinos and charginos are Wino-like. The other processes involving gauginos, 
pp W* xfxi and pp ^ Z* ^ Xi 2X1 25 are suppressed by mixing angles in the gaugino 



constrain m. This may be used to further reduce background, as the signal is peaked at m = m^^ , 
while the background is not localized in m. 

^^In the case where sleptons and fiR decay via Ir Irtti, typically Am = m^-^ — ruf-^ <C 
mj-^,mf^, and the fi is produced nearly at rest in the Ir rest frame. In the lab frame, the fi and 
Ir velocities are therefore similar, and so once sleptons cr or fiR are produced with low velocity, 
they result in a highly ionizing track, irrespective of subsequent decays. 

^^Contributions to gaugino production from t-channel squark diagrams are highly suppressed by 
large squark masses in the gauge mediation framework and are omitted in our analysis. 
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limit, and we have checked that, throughout parameter space, these are not significant 
relative to those of Eq. ([Til). The cross sections for the processes of Eq. ([TT| ) are given in 
Fig. ^ for both y/s = 1.8 TeV and 2 TeV, where here and in the rest of this subsection, we 
neglect neutralino and chargino mixings and take them to be pure Wino eigenstates. For 
equivalent gaugino and slepton masses, the gaugino cross sections are much larger, as they 
are proportional to SU(2) couplings (in contrast to right-handed slepton production, which 
is proportional to U(l) couplings in the s ^ Mz limit) and are not suppressed by the (3^ 
behavior of scalar production. 

As in Sec. |111 A[, we would like to plot cross sections for events with slow staus. The 



decays of gauginos are significantly more complicated than those of right-handed sleptons. Of 
course, the branching ratios and kinematic distributions of the cascade decays are completely 
determined once the fundamental model parameters are fixed. However, for the present 
purposes, as we will show, the cascade decays are determined to a large extent by two rather 
more meaningful parameters, namely, the physical Wino and Ir masses. In this section, 
we therefore prefer to highlight this dependence by adopting a conservative simplifying 
assumption. (The implications of these results for the fundamental parameter space will be 
discussed in Sec. |Vl| .) 

From Fig. |^, we see that for almost all parameters in the long-lived stau scenario, and 
for all parameters when > 3, m^-^ < m^o,m^±. Winos therefore decay primarily to 
left-handed sleptons, as the other possible two-body decays to Wxi, Zxi, and the "spoiler 
mode" h'^Xi are all suppressed by mixing angles in the gaugino region. The neutral Wino 
decay modes are then 

X° ^ iy'Li>' - ^ lyWhlR , (12) 

and the charged Wino decay modes are 

Xt ^ '^'I'l ^ ^%x'i - '^'I'JrIr , (13) 

where we have omitted states related by charge conjugation. From Figs. § and ^, we see also 
that the /^-X2 ^R^Xi mass splittings are never more than ~ 20% to 30% and are smaller 
than this in much of the parameter space allowed by perturbativity. In the limit that these 
states are highly degenerate, there are only two energetic final state particles in the Wino 
rest frame for any of the possible decays: the lepton from or z/^ decay, and the final state 
Ir. We therefore make the following simplification: in the rest frame of the parent gaugino, 
we assume that the only energetic particles are one lepton and the final Ir. The Ir velocity 
in this rest frame is then fixed by 

where /3 is the slow slepton's velocity and 7 = (1 — /5^)~^/^. The cross section for slow Ir 
production is then completely determined by M2 and rnj^. 

In Figs. 1^-^, we show slow Ir production cross sections for three values of mj^/M2 and 
y/s = 2 TeV, where we sum over the three slepton generations and employ the simplifying 
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assumption described above. (Again, results for = 1.8 TeV are fairly similar, and the 
corresponding figures are not presented.) To account for the roughly 27r coverage of the 
central muon upgrade detector, we have assumed an rj cut efficiency of 50%. From these 
figures, we see that for fixed gaugino mass, the cross sections are maximized for large m,- , 

''R 

that is, small /37|i.est- For M2 = 400 GeV, for example, the /37 < 0.85 efficiency is 62%, 
32%, and 7% for mj^/M2 = 0.5, 0.4, and 0.3, respectively. This can be understood as a 
consequence of the low velocity of the heavy Winos in the lab frame, which implies that the 
Ir velocity in the lab frame is mostly determined by its velocity in the Wino rest frame. 
Note that our simplifying assumption implies that we have taken the maximal possible value 
of /37|rest- In regions of the parameter space where II and X2 (and mj^ and fri^o) are non- 
degenerate at the 20% to 30% level, additional leptons can have significant energies. This 
reduces /37|i.est, and therefore increases the slow particle cross sections. Our approximation 
is therefore conservative. 

To determine the discovery reach for long-lived sleptons from gaugino production, we 
include a 75% experimental efficiency and require 5 or more events for discovery, as in 
Sec. |111 A|. The resulting discovery reaches in M2 for our representative set of luminosities 



are 



M2 = 120 to 150 GeV (v^ = 1.8 TeV, L = 110 pb"^) 

M2 = 220 to 280 GeV (v^ = 2 TeV, L = 2 fb^^) 

M2 = 270 to 340 GeV (v^ = 2 TeV, L = 10 fb"^) 

M2 = 310 to 390 GeV (v^ = 2 TeV, L = 30 fb"^) , (15) 

where the ranges are for mj^/M2 = 0.3 to 0.5. Because the signal is background-free, 
the discovery reaches grow rapidly with luminosity, and future Tevatron runs substantially 
extend our reach in SUSY parameter space. 



IV. DIMUON-LIKE EVENTS 

As discussed in the previous section, a long-lived charged slepton appears in detectors as a 
muon. Thus, even if the velocity of such a slepton is too high to be seen as a highly ionizing 
track, slepton production can lead to an excess of events with muons and an apparent 
violation of lepton universality. In this section, we discuss the simplest example of such 
an excess, namely, an apparent excess in the number of dimuon events from right-handed 
slepton pair production. 

The Drell-Yan production of long-lived fi pairs will appear as an event with essentially 
no jet activity and two high energy muon-like tracks with no ]/t- The production of right- 
handed sleptons of the first two generations may also lead to this signature. For en and 
ftR, there are three possibilities: (1) they may be stable within the detector; (2) they may 
decay via Ir — *• Irtti to leptons that are so soft that they cannot be clearly identified in 
the detector; or (3) they may decay to the NLSP and detectable leptons. In the first two 
cases, which are applicable when tan/3 is small and cr, (Ir and fi are nearly degenerate, 
the pair production of smuons and selectrons contributes to the dimuon event sample. We 
will assume this to be the case in this section. This assumption is conservative, in the sense 
that in the third scenario, cr and fiR pair production contributes to spectacular multi-lepton 
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signals, which are virtually background-free. Such signals will be discussed in the following 
section. 

The backgrounds to dimuon signals with little have been studied in the context of 



searches at CDF for new neutral gauge bosons via their decays Z' /i^/i" [Q. The 
dominant background is from genuine muon pair production pp —>■ /i+yU^.Q Before cuts, the 
cross section for this process is several orders of magnitude larger than that for slepton pair 
production. However, the signal to background ratio may be improved by requiring both 
muon-like tracks to have pt greater than some fixed p^^. In Fig. |10], we plot the cross sections 



for right-handed slepton pair and muon pair production at the Tevatron with ^/s = 2 TeV 



as a function of p^^'. Following Ref. ||2^, we require that at least one muon satisfy \ri\ < 0.6, 
and the other have |?7| < 1.0 so that its Pt may be well-measured in the central tracking 
chamber. As is evident in Fig. [1^, the pt cut is highly effective for p™^ ^ 100 GeV. 

To estimate the discovery potential for a dimuon excess, we look for an apparent violation 
of lepton universality in the ratio a{fi'^ fi~) / a{e'^ e~) . By considering this ratio instead of the 
absolute fi^fj^^ cross section, large systematic uncertainties from luminosity measurements 
and renormalization scale dependence in next-to-leading order calculations are removed. 
The statistical significance of deviations in this ratio, in units of standard deviation a, is 
given by S/ where S is the number of slepton pair events, B is the number of muon 
pair events, and we have neglected differences in e and fi acceptances and systematic errors 
in determining these acceptances. For a given slepton mass, we calculate the ratio S/ \J2B 
as a function of p™*, and optimize the choice of p™* to maximize S/\/2B. We assume an 
overall muon trigger efficiency of 75%, and sum over the three slepton generations for S. In 
Table |, we show the optimized value of S/ \/2B and the optimal choice of p'^^ for various 
values of m^^ and L. The optimal p^^ is typically p^^ ~ irif^. 

While many systematic errors are eliminated by considering the ratio a{fi^fi^)/a{e'^e~), 
systematic uncertainties in muon and electron triggering and acceptance efficiencies, of 
course, are not. The relative efficiencies may be determined, for example, by isolating 
leptons from Z decays and assuming lepton universality. However, if such efficiencies are 
uncertain at or near the level of S/B, they will seriously degrade one's ability to identify 
an excess of muons. We have therefore presented values of S/B for each m^^ in Table |; 
typical values are 10% to 20%. Relative efficiencies must be known to significantly better 
than these values for a muon excess to be detected. In the following, we will assume this to 
be the case and that the dominant uncertainty in detecting a muon excess is statistical. 

The discovery reach for a 3a excess in dimuon events is below 50 GeV for Run I, and for 
future Tevatron upgrades is 



mj^ = 70 GeV (v^ = 2 TeV, L = 2 fb"^) 
mj^ = 90 GeV (^i = 2 TeV, L = 10 fb"^) 



m, 



I = 110 GeV i^/s = 2 TeV, L = 30 fb"^) . (16) 



Of the remaining backgrounds, the largest is cosmic rays |24|. This cross section is estimated 
to be less than 1 fb for invariant dimuon masses above 200 GeV |24], and is negligible relative to 
the muon pair background if well-understood. 



12 



rrir^ 2 fb-i 10 fb^^ 30 fb^i Optimal pf^ S/B 

80 GeV L6a 3^7a 6Ict 100 GeV 22% 

100 GeV l.Ocr 2.3a 4.0c7 120 GeV 19% 

120 GeV 0.7a lAa 2.6c7 140 GeV 16% 

140 GeV 0.4a 1.0a 1.7a 160 GeV 13% 

160 GeV 0.3a 0.6a 1.1a 180 GeV 11% 

180 GeV 0.2a 0.4a 0.7a 200 GeV 9% 



TABLE I. Statistical significances of deviations in the ratio a(/i+;U^)/a(e+e^) from contribu- 
tions of long-lived charged slepton pair production to dimuon events for integrated luminosities 
L = 2, 10, and 30 fb~^ and -y/s = 2 TeV at the Tevatron. The signal from slepton pair produc- 
tion is summed over the three generations. The cut parameter pj^^ is chosen to maximize the 
significance; the optimal values used are shown. Values of S/B after cuts are as shown. 



We see that this discovery mode is clearly less powerful than the previous signal discussed 



in Sec. pT| , and we will see that it is also weaker than that of the following section. In 
fact, for an integrated luminosity of L = 2 fb~^, it is not possible to achieve a 3a excess, 
given current bounds from LEP. However, this mode may still be useful as a supplementary 
signal. For example, if highly ionizing tracks consistent with a certain charged particle mass 
are seen, the discovery of even a small but quantitatively consistent dimuon excess could be 
an interesting confirmation of the supersymmetric interpretation of these exotic signatures. 
Note also that it may be possible to improve the sensitivity of the dimuon channel by 
loosening the t] cuts, otherwise improving acceptance, or by combining data from the two 
detectors. 

Given the rather weak results for the Tevatron, we consider the possibility of observing a 
dimuon excess from pp —>■ IrI'^ at the LHC with ^/s = 14 TeV. We require r] < 2.2 for both 



sleptons, corresponding approximately to the muon coverage of both LHC experiments [25 



and, as above, assume a 75% experimental detection efficiency and sum over 3 generations. 
In Table m we present results for integrated luminosities of 10 fb~^ (one year, low luminosity), 
100 fb~^ (one year, high luminosity), and 600 fb~^ (corresponding to multi-year runs and/or 
combined data sets from the two detectors). Not surprisingly, the results are much improved. 
For example, for L = 100 fb~^, a 3a excess will be discovered for rrij^ ^ 300 GeV, where 
we have again assumed that systematic uncertainties in acceptance are small compared to 
statistical uncertainties. Of course, at the LHC, a number of other signals of gauge-mediated 
supersymmetry are likely to be evident, and even squark and gluino cross sections may be 
substantial. 



V. MULTI-LEPTON SIGNALS 

If stau NLSPs are misidentified as muons, supersymmetric particle production can lead 
to events that appear to have large numbers of leptons in the final state. Here we focus 
on events with 5 or more isolated leptons and very little hadronic activity. Assuming that 
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ruf^ 10 fb"^ 100 fb"^ 600 fb"^ Optimal p^^^ S/B 

100 GeV Slo 26a 63a 150 GeV 28% 

200 GeV 2.4a 7.7a I9a 250 GeV 20% 

300 GeV l.lo- 3.4a 8.3a 350 GeV 16% 

400 GeV 0.5a 1.7a 4.1a 450 GeV 13% 

500 GeV 0.3a 0.9a 2.3a 550 GeV 11% 

600 GeV 0.2a 0.5a 1.3a 650 GeV 9% 



TABLE II. Same as in Table |, but for integrated luminosities L = 10, 100, and 600 fb ^ and 
^ = 14 TeV at the LHC. 



the probability for jets and photons to fake leptons is well under control, the backgrounds 
to such events are small, and so even a few signal events may constitute a discovery signal. 
Such multi-lepton states may result from either right-handed slepton or gaugino production, 
and, as in Sec. |T|, we consider these in turn. 



A. Slepton Production 

For small tan (3, right-handed slepton pair production leads to dimuon events, as analyzed 
in Sec. However, for large tan/3, the intergenerational splitting from the tau Yukawa 
coupling may be substantial, and selectrons and smuons may decay to NLSPs and sufficiently 
energetic leptons to be detected as multi-lepton events. 

To understand what values of tan/3 are required for this scenario, let us consider the 
various contributions to the intergenerational mass difference. Left-right stau mixing gives 
an important contribution, 

In addition, the splitting receives contributions from RG evolution. Roughly speaking, this 
effect is proportional to ln(M/MMssM) and the beta function for M?^, and hence may be 
estimated as 

^^^^ ^ 2^ ^ ^ ^ ["^^^ ln(M/MMssM) 



X O 



4 

Tfl 

ln(M/MMssM: 



X AmLR , (18) 



where Ur is the r Yukawa coupling constant, v = 246 GeV, and in the last line we have 
assumed large tan p. AmRc niay be sizeable compared to Attt-lr when nif^ is large relative 
to the electroweak scale. In any case, both contributions reduce m^^ and they are both 
enhanced when tan/5 is large. As a result, mass differences between fi and the other right- 



handed sleptons can become large for high tan/?. In Fig. 11, we plot the mass splitting 
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Am = rriej^ — rrif-^ ~ m^^ — m^^ as a function of m^^ for fixed = 3, M = 10^ GeV, /i > 0, 
and various tan/3. If tan/3 is small, Am < rrir, and all three right-handed sleptons are stable 
in the detector. However, for tan /3 = 30, for example, mass splittings of 40 GeV or larger 
are possible, and and (Ir may decay inside the detector to electrons and muons that are 
likely to be observable. 

We now consider the scenario in which right-handed selectron and smuon pair production 
leads to final states (ee, /i/i)rrfifi. For high tan/3 and correspondingly large mass splittings 
Am, the resulting electrons and muons, and even the leptonic decay products of the tau 
leptons, are typically energetic enough to be detected. Demanding that at least one r decay 
leptonically, this leads to signals of either 5 energetic, isolated leptons and the hadronic 
decay products of a r, or 6 energetic, isolated leptons with no significant hadronic activity. 
We expect the backgrounds to such signals to be extremely small. On the other hand, for 
more moderate values of tan/3, the leptons are typically softer, and the number of events 
with 5 or 6 observable leptons is reduced. While it may be possible to lower the requirement 
on lepton energy, significant backgrounds may then enter. 

To obtain an estimate of the possible reach of this mode, we assume, for simplicity, that 
tan P is sufficiently large, for example, tan /3 ^ 30, such that to a good approximation, all 
leptonic decay products are energetic enough to be detected. We take the signal to be 5 or 
more isolated leptons with no significant hadronic activity, other than that consistent with 
the hadronic decay products of a r. Slepton pair production then contributes to the signal 
when at least one r decays leptonically. With these assumptions, we expect the signal to 
be essentially background-free and we therefore require a discovery signal of 5 events. To 
roughly incorporate the effects of detector acceptances, we include an overall experimental 
efficiency of 50%. From Fig. ^, but for 2 generations, we find the following discovery reaches 
for slepton masses: 

m^^ = 60 GeV (^i = 1.8 TeV, L = 110 pb"^) 
mj^ = 130 GeV (v^ = 2 TeV, L = 2 fb"^) 
mj^ = 170 GeV (v^ = 2 TeV, L = 10 fb-^) 

mj^ = 210 GeV (v^ = 2 TeV, L = 30 fb"^) . (19) 

As noted above, these estimates of the discovery reach are valid only for large tan /3, 
and are degraded for lower tan/5. A quantitative determination of the dependence of the 
reach on tan requiring a careful study of backgrounds and the optimization of cuts, merits 
further study, but is beyond the scope of the present analysis. 



B. Gaugino Production 



As discussed in Sec. [Ill B| , in the stau NLSP scenario, Wino-like charginos and neutralinos 
typically decay through left-handed sleptons, with cascade decays given in Eqs. ([T2|) and 
(p!3|). Assuming that right-handed sleptons appear in the detector as muons, the final state 
therefore contains many charged leptons and no jets. For scenarios with gauginos that are 
highly degenerate with sleptons, some of the real leptons in the final state may be too soft 
to detect. However, this applies only to points extremely close to the m^^ = Mi boundary 
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of Fig. 0, and we will assume that this is not the case. Gaugino production then results in 
the spectacular signal of from 5 to 7 charged leptons.0 

As in Sec. |V A\ , we take the signal to be 5 or more isolated leptons with little hadronic 
activity. We demand 5 events for a discovery signal, and again include an overall exper- 
imental detector efficiency of 50%. For the 7 lepton events, this is likely to be a rather 
conservative estimate, as only 5 or more need be detected. We may then determine the 
discovery reach for multi-lepton signals in this scenario from Fig. ^. The discovery reaches 
for our representative data samples are 



M2 


= 190 GeV 


(v^ = 


1.8 TeV, L 


= 110 pb 


M2 


= 310 GeV 


(v^ = 


2 TeV, L = 


2 fb-i) 


M2 


= 370 GeV 


(v^ = 


2 TeV, L = 


10 fb"^) 


M2 


= 410 GeV 


(v^ = 


2 TeV, L = 


30 fb"^) 



The multi-lepton signal from Wino production is thus one of the most sensitive probes, 
and probes chargino and neutralino masses far into the ranges typically expected for Wino 
masses. 



VI. COMPARISON OF RESULTS 

In the previous 3 sections, we have considered a variety of SUSY discovery signals in the 
long-lived stau NLSP scenario. They are: 



(A) Highly ionizing tracks from slepton pair production (Sec. |III A| ) 



(B) Highly ionizing tracks from Wino pair production (Sec. |IIIB|) . 

(C) An excess of dimuon events from slepton pair production (Sec. [TVI ). 

(D) Multi-lepton signals from slepton pair production (Sec. [V A|) . 

(E) Multi-lepton signals from Wino pair production (Sec. |VB| ). 

For each of these signals, we have obtained discovery reaches in terms of physical right- 
handed slepton and/or Wino masses by the analyses described above. We have found that 
the discovery reach for slepton masses from (A), (C), and (D) are numerically weaker than 
those for gaugino masses from (B) and (E). However, in typical gauge mediation scenarios, 
the Winos are predicted to be substantially heavier than the right-handed sleptons. It is 
therefore of interest to determine what the relative strengths of these probes are in terms of 
the fundamental parameter space of the gauge mediation model. 

In Figs. O-ITB, we summarize all of our previous results by presenting discovery reaches 



in the fundamental parameter space for various integrated luminosities. We display these 



^^We assume here that e/j and /i/j do not further decay through Iji Irtti. If such decays result 
in detectable leptons, the signals are even more spectacular. 
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reaches in the (M, A) plane, and fix the remaining three parameters tan (3 = 3, /i > 0, and 



N^. The number of messenger representations is A^s = 3 for Figs. |12| and |14|; in Figs. |13| and 
|T5|, we illustrate the dependence of these results by taking = 4. For each point in 
the plane, the complete superpartner spectrum is specified, and, in these figures, D-terms 
for scalar masses and all mixing effects are included in calculating physical masses and 
production cross sections. In each figure, we shade the area in which > m^o; in the 
remaining region, the stau is the NLSP. 

To guide the eye in relating the fundamental model parameters to more physical param- 
eters, we plot dotted contours of constant m^^. Also, since Mi/gf is invariant under one- loop 
RG evolution and the weak scale coupling constants 5'i(MMssM) are known, the weak scale 
gaugino masses Mj(MMssM) are proportional to A, independent of M. In particular, the 
SU(2)i gaugino mass is given by 

M,(MMssM) = Ar,A^M^^. (21) 

We therefore also give vertical axis labels for M2(Mmssm) with Mmssm = 1 TeV. This axis 
gives approximately the physical masses of the Wino-like chargino and neutralino. 

For each point in parameter space with a stau NLSP, we determine the potential for 
SUSY discovery for each signal according to the analyses described in Sees. PTTHVI. Contours 
for signals (A), (C), and (D) therefore roughly follow contours of constant right-handed 
slepton mass, and contours for signal (E) roughly follow contours of constant chargino and 
neutralino mass. For signal (B), the contours also roughly follow chargino and neutralino 
mass contours, but rise slightly relative to these contours as M increases: in this direction, 
m,- increases, and the signal is strengthened, as discussed in Sec. [IIIB. 



We begin by presenting results in Figs. O and IT3| for the present Tevatron data sample 



of L = 110 pb^^ at ^/s = 1.8 TeV. Several comments are in order. First, we have assumed 
tan/5 = 3, and so slepton production does not result in multi-lepton events. Signal (D) is 
therefore not applicable. Second, signals (A) and (C) are too weak to appear in the plots. 
These signals are both from slepton production, and are essentially dependent on only the 
right-handed slepton mass. Given only the current Run I data, they are therefore weaker 
than the current LEP bound of m^^ > 75 GeV, independent of model assumptions. For 
signal (B), the reaches plotted are conservative, as discussed in Sec. |111 B| . Finally, signal 
(E), rnulti-lepton events from Wino production, requires that the leptons from W II and 
B ^ Ir transitions be energetic enough to be detected. For A^s = 3 and 4, this is always 
true for the former decay in the stau NLSP region. However, the latter decay may result in 
soft leptons, and the bound from signal (E) is therefore weakened for points very close to 
the boundary of the shaded region. We have not included this effect in the figures. 

For the particular parameters we have chosen, a number of interesting conclusions may 
be drawn. Given the current Run I data, the slepton production signals (A) and (C) are 
superseded by the requirement that the parameters yield the stau NLSP scenario in the 
first place. Thus, for low tan/3, only the gaugino production signals (B) and (E) provide 
non-trivial probes of the parameter space. However, even given only the current Tevatron 
data sample, signals (B) and (E) do probe regions of parameter space not excluded by the 
current LEP bound nif^ > 75 GeV |T^. For = 3, we see from Fig. |T2| that signal (B) 
is sensitive to parameter regions beyond this bound, and signal (E) probes regions with 
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^ 110 GeV, which are even beyond the ultimate kinematic reach of LEP II. For = 4, 
the ratio mfj/M2 is reduced, and the strength of gaugino production mechanisms relative 
to slepton mechanisms is slightly weakened. However, as seen in Fig. |T3|, signal (B) is still 
competitive with the current LEP bound, and signal (E) will significantly extend the reach 
of LEP II for certain values of M. 

In the context of the stau NLSP scenario, if no signals are found in analyses of the 
available Tevatron data, current bounds on Wino masses will be extended to roughly 200 
GeV. In addition, in the framework of the particular model we have described, an indirect 
bound on m^^ may be set. For example, for = 3, we see from Fig. that stau masses 
below ~ 90 GeV would be excluded by the absence of signal (E). Of course, such indirect 
bounds on stau masses rely heavily on the specific model assumptions described in Sec. fl[ 

In Figs. |I4| and [1^, we present results for the same model parameters, but for future 
Tevatron runs with a/s = 2 TeV and L = 2, 10, and 30 fb~^. Again, as we have assumed 
tan/3 = 3, discovery reaches from (D) multi-lepton signatures of slepton pair production are 
inapplicable. Signal (C) is again weak relative to the others. For integrated luminosities 
of 2 fb~^ or more, direct probes of slepton masses from signal (A) and indirect probes 
from signals (B) and (E) all exceed m^^ ~ 100 GeV, and SUSY discovery in these modes 
is possible, even given potential LEP II bounds on the stau mass. For L = 2 fb-^ the 
multi-lepton gaugino production signal (E) is dominant, and indirectly excludes stau masses 
below ~ 120 GeV. However, for increasing luminosity, the relative importance of the slow 
slepton signal (A) increases. For L = 10 fb~^ and L = 30 fb~^, signal (A) becomes the most 
sensitive probe, with sensitivity to regions of parameter space with Wino masses in excess 
of 400 to 600 GeV. 



VII. CONCLUSIONS 

Gauge-mediated supersymmetry breaking models provide an elegant solution to one of 
the most puzzling issues in supersymmetry, namely, the supersymmetric flavor problem. 
The phenomenology of these models is to a large extent dominated by the character of the 
NLSP. In the most popular gauge mediation scenarios, there are four options, as the NLSP 
may be either the neutralino or stau, and may be either stable or unstable in detectors. 
Three of these scenarios have received fairly extensive attention in the past. In this study, 
we have given for the first time a detailed treatment of prospects for SUSY discovery at the 
Tevatron in the fourth scenario, in which the NLSP is a long-lived stau. As discussed in 
Sec. H, this scenario is realized in a wide range of parameter space in models with messenger 
sectors with iVs > 2. 

In such scenarios, all superparticles decay ultimately to the lighter stau fi, which ap- 
pears in collider detectors as a (possibly slow) muon. We have analyzed the prospects for 
discovering SUSY by detecting highly ionizing tracks from slow fi production, an apparent 
excess in dimuon events from stau pairs, and the appearance of spectacular multi-lepton 
events from superpartner pair production followed by decays. 

Of the signals considered, the dimuon excess is least promising, as it suffers from the 
large background of genuine /i"*"/^" production. Its use appears to be limited to confirming 
the supersymmetric interpretation of other, more sensitive signals. Note, however, that a 
dimuon excess is only the simplest example of an apparent violation of lepton universality 
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from long-lived fi production. Many other possibilities exist. If a statistically significant 
excess of /i over e events appears in any event sample, we strongly encourage the investigation 
of the long-lived charged slepton hypothesis, either by looking for highly ionizing slow tracks 
in this event sample, or by forming distributions of kinematic variables under the assumption 
that the muons are in fact sleptons of some given mass. 

The remaining signals are sufficiently spectacular that they are essentially background- 
free after the cuts we have detailed in Sees. ^ and |V|. We have analyzed various data sam- 
ples, beginning with Run I data with ^/s = 1.8 TeV and integrated luminosity L = 110 pb~^. 
With this current data, we find that the sensitivity from slepton production is not compet- 
itive with the current bound on long-lived staus of m^^ > 75 GeV from combined LEP 
analyses. However, Winos with masses beyond the LEP kinematic limit may be produced 
in significant numbers at the Tevatron. Wino pair production therefore probes regions of 
parameter space beyond the current LEP bound, and even beyond the possible ultimate 
kinematic reach of LEP II. The analysis of currently available Tevatron data therefore pro- 
vides interesting new probes of the gauge mediation parameter space. 

Given the great interest in future Tevatron upgrades, we then considered what improve- 
ments may be expected from future Tevatron data. As the backgrounds to the signals we 
consider are extremely suppressed, we find that significant improvements may be expected. 
For future Tevatron runs with integrated luminosities of 2, 10, and 30 fb~^, we estimate 
that SUSY may be discovered for right-handed slepton masses of 110, 180, and 230 GeV 
and Wino masses of 310, 370, and 420 GeV, respectively. These new and spectacular sig- 
nals are therefore sensitive to much of the typically expected superpartner mass range and 
provide a promising avenue for SUSY searches at Tevatron runs in the near future. 
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M (GeV) 

FIG. 1. Contours of Mf^(MMSSM) = Mi(Mmssm) for Mmssm = 1 TeV and tan/3 = 3, 10, and 
30. In the region above the contours, fi is the NLSP; in the region below, Xi is the NLSP. In 
the shaded region, gauge couphng constants become non-perturbative below the GUT scale under 
two-loop RG evolution. 
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FIG. 2. Contours of the ratio of soft SUSY breaking parameters Mf^(MMSSM)/-^i(-^MSSM) 
for tan/3 = 3 and A/mssm = 1 TeV. In the lower shaded region, Aff^ (Mmssm) > -^i(-Wmssm)- In 
the upper shaded region, gauge coupling constants become non-perturbative below the GUT scale. 
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FIG. 3. Contours of the ratio of soft SUSY breaking parameters Mi^{Mmssm)/M2{Mmssm) for 
tan/3 = 3 and Mmssm = 1 TeV. In the lower shaded region, Mf^(MMSSM) > -^iCMmssm)- In the 
upper shaded region, gauge couphng constants become non-perturbative below the GUT scale. 
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FIG. 4. Cross sections for pp ^ j* , Z* fiff for ^/s = 1.8 TeV (lower) and 2 TeV (upper), 
and fi ffl. 
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FIG. 5. Cross sections for the production of at least one slow fi from — ^ 7*, Z* fif* at 
y/s = 2 TcV, where we have required |r/| < 0.6 for the slow fi. Contours correspond to P'j < 0.4 
(dotted), 0.7 (dashed), 0.85 (solid), and 00 (dot-dashed). 
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FIG. 6. Cross sections for pp W* xfx2 (solid) and pp 7*,Z* xfxf (dotted), for 
^/s = 1.8 TeV (lower) and 2 TeV (upper). Chargino and neutralino mixing effects are neglected. 
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FIG. 7. Cross sections for slow Z/j production from chargino and neutralino production with 
m-^^/Al2 = 0.5, and summed over the three slepton generations. Contours correspond to (3^ < 0.4 
(dotted), 0.7 (dashed), 0.85 (sohd), and oo (dot-dashed). An efficiency of 50% is included for the 
requirement that the slow slepton have |?7| < 0.6. 
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FIG. 8. Same as in Fig. 0, but for /M2 = 0.4. 
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FIG. 9. Same as in Fig. |^ but for mj /M2 = 0.3. 




(GeV) 

FIG. 10. The cross section for fi pair production at the Tevatron after requiring pT > Pj^^ for 
each stau, as a function of pj^^. We require one fi to have |r/| < 0.6 and the other to have \ri\ < 1.0, 
and we take ^/s = 2 TeV and ruf^ = 80, 100, 120, 140, 160, and 180 GeV (sohd lines from above). 
The cross section for the background pp — > /^"'"/i" is also shown (dotted). 
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200 



FIG. 11. The mass splitting Am = mg^ — nif^ ~ mp_j^ — mf^ as a function of mfj for fixed 
= 3, M = 10^ GeV, /x > 0, and the tan P values shown. 
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FIG. 12. Summary plot for = 1.8 TeV and L = 110 pb"^ (Run I), and fixed N5 = 3, 
tan/3 = 3, and /j, > 0. Solid contours give the discovery reach for (B) highly ionizing tracks from 
gaugino production, and (E) multi-lepton signals from gaugino production. Contours of constant 
ruf^ are given by the dotted curves. In the shaded region, fi is not the NLSP. 
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FIG. 13. Same as in Fig. 12, but for iVs = 4 



28 



600 



500 - 



> 

S 400 



300 



200 



500 



2, 400 



300 



200 



500 - 



> 



400 - 



300 



200 



rTTTTTir rTTTTTnl — 'L 

1 1 1 Mill 1 1 1 Mill r 1 1 Mill 

: (a) ' • 


TTTTTTir rTTTTl 
1 1 llllll 1 1 1 III 

2S0'"-. 


TTTTTTTm — T 
* J 1 1 Mill 1 

300'-. 


1 1 Mill 1,1 1 Mill 

350 ^ 


200"-.. 






'•.. - 


ISO'--. 








E 








_ 1 1 .1 1 


B 7 
mwi Mini 


1 1 


1 


:..(b) ■ 


2S0''-, 


300'-. 


350 _i 


E 








ISO'--.. 








" C ' °° 

^MT™hTTtttlli li 


'■ ^1 


1 1 
1 1 


1 

llllllll .1 llllll 


- (c)\ • 

r • . A\ 


250''-, 


300'-. 


350 _i 


200"-.. 

::_ E 














ISO'--.. 


B 1 














J---....100 

1 1 iiiiiil 









- 70 



^ 60 



H 50 H 

CD 
< 



— 40 



30 



- 70 



^ 60 



^ 50 



— 40 



■-^ 30 



70 



^ 60 



^ 50 



— 40 



■-^ 30 



10' 10= 10' loMo'lO '10 ^^0 Mo 
M (GeV) 



CD 
< 



CD 
< 



FIG. 14. Summary plot for -y/i = 2 TcV and integrated luminosity (a) 2 fb^^, (b) 10 fb^-*^, and 
(c) 30 fb~^. We fix = 3, tan/? = 3, and > 0. Solid contours give the discovery reach for 
(A) highly ionizing tracks from slepton production, (B) highly ionizing tracks from gaugino pro- 
duction, (C) a dimuon excess from slepton production, and (E) multi-lepton signals from gaugino 
production. Five events are required for (A), (B), and (E), and a 3a excess is required for (C). 
Contours of constant mf^ are given by the dotted curves. In the shaded region, fi is not the NLSP. 



29 




30 



